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Amendments to the Specification 

Please replace the section beginning at page 3, line 12, with the following 
rewritten section: 

BRIEF DESCRIPTION OF THE DRAWINGS 
For a more complete understanding of the present invention and for further 

features and advantages, reference is now made to the following description, taken in 

conjunction with the accompanying drawings, in which: 

FIGURE 1 is a flowchart illustrating an example method for generating a 

microarray; 

FIGURE 2 illustrates background noises and includes FIGURES 2(a), 2(b), and 

2(c) ; 

FIGURE 3 illustrates noise settings and includes FIGURES 3(a) and 3(b) ; 

FIGURE 4 illustrates an example cDNA microarray spot model; 

FIGURE 5 illustrates variability in spot size and spread from the spot size and 
includes FIGURES 5(a). 5(b). and 5(c) ; 

FIGURE 6 illustrates inter-spot grid spacing and includes FIGURES 6(a). 6(b), 
and 6(c) ; 

FIGURE 7 illustrates an effect of a radius drift variation and includes FIGURES 
7(a). 7(b). and 7(c) ; 

FIGURE 8 illustrates chord rate settings and includes FIGURES 8(a). 8(b), and 

8{c}; 

FIGURE 9 illustrates an edge noise of spots and includes FIGURES 9(a), 9(b), 
and 9(c) ; 

FIGURE 10 illustrates fluorescent detection response characteristic functions and 
includes FIGURES 10(a). 10(b). 10(c), and 10(d) ; 

FIGURE 1 1 illustrates possible scatter plots due to various response conversions 
for different fluorescent channels and includes FIGURES 1 1(a). 1 Kb), and 1 1(c) ; 
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FIGURE 12 illustrates increased spike noise levels and includes FIGURES 12(a), 

12('b).and 12(c) ; 

FIGURE 13 illustrates scratch noise and parameter settings and includes 
FIGURES 13(a). 13(b). and 13(c) : 

FIGURE 14 illustrates parameter settings for snake noise and includes FIGURES 
14(a). 14(b). and 14(c) ; 

FIGURE 15 illustrates convolution kernels and includes FIGURES 15(a) and 

15(b) ; 

FIGURE 16 illustrates a three-dimensional profile before and after smoothing mid 
includes FIGURES 16(a). 16(b). and 16(c) ; 

FIGURE 17 illustrates a full size array simulation with different parameter settings 
and includes FIGURES 17(a) and 17(b) ; 

FIGURE 18 illustrates a comparison between a simulated signal versus an 
extracted signal from a microarray image analysis program and includes FIGURES 18(a) 
and 18(b) ; and 

FIGURE 19 illustrates simulated images exhibiting undesirable noise conditions 
and includes FIGURES 19(a). 19(b). 19(c). and 19(d) . 

Please replace the paragraph beginning at page 1, line 14, with the following 
rewritten paragraph: 

Since the inception of cDNA microarray technology [1] as a high throughput 

method to gain information about gene functions and characteristics of biological 
samples, many applications of the technology have been reported [2-10]. With the 
improvement of the technology, including fabrication, fluorescent labeling, hybridization, 
and detection, many computer software packages for extracting signals arising from 
tagged mRNA hybridized to arrayed cDNA locations have been designed and applied in 
various experiments [11 13] [11-12] . As reported in [1 1], a target detection procedure has 
been implemented that utilizes manually specified target arrays, extracts the background 
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via the image histogram, predicts target shape by mathematical morphology, and then 
evaluates the intensities from each cDNA location and its corresponding ratio quantity. 



Please replace the paragraph beginning at page 7, line 27, with the following 
rewritten paragraph: 

Rather than be constant across the entire microarray, the mean of the background 
noise may vary owing to various scanning effects. It can take different shapes: parabolic, 
positive slope, or negative slope. In this case a function g(x, y) is first generated 
(parabolic, positive slope, or negative slope) to form a background surface and normal 
noise is added to it pixel wise. Thus, the background intensity is of the form h ~ N(p,b, 
0/) with \ib = jg(x, y), where y ~ Uia, b) is the targeted background noise level. 
Background deviation is set independently for each channel: = k^^Ph andcj^^ = '"^ 
FIGURE 2 shows various noiac backgrounds background noises with the background 
deviation factor set at = k,^ = 0.1 00%) . The sienal-to-noise ratio (SNR) is s et at 1.0 
for the slides of FIGURE 2. The slides have the following set tings: FIGURE 2(a), 
Parabolic background noise. FIGURE 2(b). Positive slope background, FIGURE 2(c), 
Negative slope background all with global noi se parameter. 



Please replace the paragraph beginning at page 8, line 4, with the following 
rewritten paragraph: 

In many practical examples, the non-specific hybridization at the target location 
may be different from its peripheral region. Although one may have trouble pin-pointing 
this particular observation under normal conditions owing to signal interference, it is 
sometimes unmistakable when locations assumed to be weakly expressed, or not 
expressed at all, carry some non-zero readouts, or the intensity in the center is stronger 
than the doughnut-ring if the printed target is doughnut-shaped. This artifact is simulated 
under a gradient noise condition by allowing the background for the center holes to be at 
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higher levels than the signal intensities. Hence, there is an option to use global 
background or local background information to set the noise parameter for the center 
hole. FIGURE 3 shows the effects of using local and global background parameters. Ihe 
examples of FIGURE 3 show different noise settings for a spot's inner hole. FIGURE 
3fa) uses a global background parameter to fill the center hole. FIGURE 3(b) uses local 
background for filling the center hole. The background noise is set to sloped type with 
SNRof 1.5. This effect may not appear everywhere in a simulated image; however, it is 
often sufficient to require appropriate procedure design in the image analysis program to 
lessen the penalty. The effects of weak targets will be fiirther studied in later sections. 

Please replace the paragraph beginning at page 9, line 8, with the following 
rewritten paragraph: 

Prior to distortion and noise, the cDNA deposition spot is considered to be circular 
with random radius S. The mean of the radius is set according to the array density and its 
variance relates to the consistency of spot size. S is modeled by a normal distribution 
having mean [is and variance a/, S ~ Ni\is, a,), with the standard deviation being a pre- 
determined proportion, k^, of the mean, or 5 ~ NiHs, k[is). The radius mean is set for 
every block, and randomized over a small range within the array. The block randomness 
of \is is modeled by a uniform distribution, ~ Uis^, Sb). FIGURE 5 shows parts of 
blocks with spot radii depending on the number of spots in a block. For FIGURES 5a 
through 5c, the block portions aro for block oizos (10, 15), (25, 15), and (25, 15), 
ror,poctivelv FIGURES 5(a). 5(b), and 5(c) show the variability in spot size and spread 
from its size. The spot radius distribution is automatically set depend ing on the number 
of spots in a block (width, height), where (col, row) denotes the number of columns and 
rows within the block, respectively. In FIGURE 5(a). the block portion is for size 
(10.15). u.~Uf23.3 24.31. In FIGURE 5(b). the block portion is for size (20, 25) , 
~U[12.6 13.6]. In FIGURE 5(c). the block portion is for size (25. 45) . u.~Ur5.45 
6.451. Standard deviation L. equals 1%. 7%. and 20% of radius, respectively. 
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Occasionally, a spot overlaps with it neighbors (FIGURE 5o) (FIGURE 5(c)) when ks 
is set to a larger proportion. This situation simulates the condition where too much cDNA 
solution is deposited and/or the drying process may be slow in comparison to the liquid 
spreading process. 

Please replace the paragraph beginning at page 9, line 22, with the following 
rewritten paragraph: 

Depending on the robot arm and printing ability of the pins, the inter-spot 
distance, Gsp, may vary. Owing to the physical mechanics of the robot arm, the block size 
(pixel units) is fixed in most cases. The inter-spot distance can be set to accommodate 
spot size and random variation in spot radii. The effects are illustrated in FIGURE 6, 
where the number of rows and columns are fixed. FIGURE 6 shows inter-spot grid 
spacing. FIGURE 6(a) shows G.^. = 3 pixels. ^u~Ur9.5 10.51. FIGURE 6(b) shows 
G^ . - 6 pixels, u ^ ~Ur8 91. FIGURE 6(c) shows G^. = 10 pixels, u^~U[6.5 7.51. The 
example shown has (35. 20) rows, columns respectively with k, = 0.05. 

Please replace the paragraph beginning at page 10, line 8, with the following 
rewritten paragraph: 

Some microarray scanners capture two fluorescent signals in two passes of 
scarming. Due to the mechanical homing error, the two fluorescent channels may not 
align exactly. In these settings, some small offset between the two channels can be 
observed. This offset may occur at sub-pixel resolution. To simulate this offset, the model 
offers a random offset between the centers of the two channels. It is achieved by 
randomly offsetting the spot center of the second chaimel by one pixel in either of the 
horizontal and vertical directions. These offsets are applied following application of the 
spot drifts in the first channel. FIGURE 7 illustrates the spot drift. FIGURE 7 shows the 
effect of radius drift variation (Pn , Rji, Rai ). The settings for FIGURE 7(a) are (0.05, 5, 
100). The settings for FIGURE 7(b) are (0.25. 15. 100). The settings for FIGURE 7(c) 
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are ro. 5. 50. 100 V As the activation probability with drift ranee is set higher, a spot drifts 
away from its center. 

Please replace the paragraph beginning at page 12, line 1, with the following 
rewritten paragraph: 

Once the number of chords for a spot is determined, the distance, L, of each chord 
center to the edge is selected from a beta distribution: L ~B(ai, pi). Inter-block 
variability is modeled by allowing and pi to be randomly selected from uniform 
distributions: ~ U{a^, ba), and ~ U{a^, b^). Owing to the large family of shapes 
generated by beta distributions, this provides a wide range of distributions for L. Finally, 
the chord locations are chosen uniformly randomly according to an angle 9 ~ U(0, 2n). 
FIGURE 8 illustrates the effect of selecting increased chord rates. FIGURE 8 shows 
different chord rate settings for each of the slides. T he probability weights for (0,1,2,3,4) 
chords were set at the following levels : FIGURE 8 (a) = 0.70,/?, = 0.30;£2.^0i^3^ 
p ^ = 0- FIGURE 8 (b) po = 0.20, p, = 0.40, p2 = 0.25, = 0. 1 5 ; FIGURE 8(c) po = 

0,px = 0.l0,p2 = 0.40, p3 = 0.30, p4 = 0.20. The chord rate is reset at the beginning of a 
block. 

Please replace the paragraph beginning at page 13, line 9, with the following 
rewritten paragraph: 

where 5 controls the threshold and hence the edge noise, and A denotes the symmetric 
difference. 8 is used as controlling parameter. S' is a binary mask giving the spatial 
domain of the spot. FIGURE 9 shows edge noise for various 6 thresholds. FIGURE 9 
shows the edge noise on the spots. A noise controll ing parameter (5) can be set from [0, 
1.01. The examples of FIGURE 9 show an increase d edge noise effect. For FIGURE 
9fa). 5 - 0.25. For FIGURE 9(h\ 5 = 0. 1 . For FIGURE 9(c), 5 = 0.03. 
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Please replace the paragraph beginning at page 13, line 20, with the following 
rewritten paragraph: 

It is well known that the distribution of gene expression levels within a cell closely 
follows an exponential distribution ^ £13]. Given a microarray containing N genes, the 
intensity levels , for A; = 1 , . . ., assumed to be related to the expression levels of iV 
genes, are simulated by an exponential distribution. This intensity level Ik is considered to 
be the ground-truth signal that is not directly measurable from microarray, since from 
either biological or bio-chemical processes, from mRNA extraction up to the 
hybridization process, some variation will be introduced into measurement of final 
fluorescent signal strength. For each microarray, a particular exponential distribution with 
mean p is first chosen (for a detection system with gray-level up to 65,535, p is usually 
selected around 3000). Then at each spot location, which is assumed to represent one 
unique gene, one ground-truth signal level h is generated from the exponential 
distribufion. For two observable measurements (Rk, Gk) from two fluorescent channels, 
two numbers are generated from a normal distribution with mean of 4 and standard 
deviation of a/fc, where a is a pre-determined coefficient of variation, which is usually 
about 5% to 30% depending on the assumed biological relation between the two 
channels. 

Please replace the paragraph beginning at page 15, line 3, with the following 
rewritten paragraph: 

Owing to various reasons, such as imprecise quantities of starting mRNA for the 
two channels, different labeling efficiencies, or uneven laser powers at the scanning 
stage, in actual microarray experiments there may not be equal intensities even if two 
channels use exactly the same labeled mRNA. Moreover, one may not be able to assume 
that the fluorescent intensity is linearly related to the expression level. In fact, it is very 
difficult to determine the exact form of the response function from expression level to 
intensity due to the complex combination of bio-chemistry to photon-electronics. A 
family of fimctions that covers most of the understandable conditions, shown in FIGURE 
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10, such as delayed response, saturation (which is an embedded feature in the digital 
system since in general no gray-level can pass 16-bit binary digits in a typical microarray 
system), and unbalanced channel intensity, is selected. FIGURE 10 shows fluorescent 
detection response characteristic functions. In the figures, the middle (circled) curve is 
the reference fimction with the parameters of (ao, ai, a?, Qy,) ^ (0, 100, -1, !)• Also, in all 
the figures, the x-axis is the input signal intensity, and v-axis is the observed signal 
intensity, and both are in log in -scale. FIGURE 10(a) sh ows delayed response at various 
levels, with fixed a^j = Q and a. = 1. FIGURE 10(b) shows different amplification levels, 
with fixed a n = 0 and a9 - -1 ■ FIGURE lQ(c) shows different response curvature, with 
fixed an = 0 and a. = 1 ■ FIGURE 10(d) shows some other parameter settings, with fixed 
a, = 1. This simulation is intended to facilitate understanding as to what is the best way 
for expression ratio normalization, whether linear based methods will be sufficient or 
non-linear based methods may be needed. The function family is characterized by four 
parameters, (ao, ai, «2, ay), and the function form is given by 

Please replace the paragraph beginning at page 16, line 15, with the following 
rewritten paragraph: 

The scatter plots in FIGURE 1 1 show the effects of the channel normalization. By 
choosing different parameter sets, one can simulate many of the situations observed in 
real microarray images. FIGURE 1 1 shows possible scatter plot due to various resp onse 
conversions for different fluorescent channels. 10.000 data points (gene e xpression 
levels) were generated bv the exponential distribution with m ean of 3000. After passing 
through two fluorescent channels (with some respo nse characteristic fiinctions as shown 
in FIGURES 1 1(a) to 1 1(c)). data variations were added by passing each data point 
through a normal distribution with the standard deviation to b e 15% of mean expression 
signal. FIGURE 1 1(a) shows without anv alteration (or equi valentlv. set parameters for 
the response function to be (an . a., a,, a, ^= (0.1. -1. U). and assume t h e signal intensities 
from red channel and green channel are equivalent (a simula ted self-self experiment). 
FIGURE ] 1 (b) shows a banana shape. Intensity in the green channel pass a response 
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function with parameters (a n. a u ao- a^>= (0. 500. -1.0. where the red channel takes the 
parameters (0, 10, -1. 1). FIGURE 1 1(c) shows a sinusoid shape. The red channel's 
response function with parameters (0. 100'^° \ -0.7. I) and the green channel with (0. 
lOO'^^l -0.9. 1). 

Please replace the paragraph beginning at page 17, line 7, with the following 
rewritten paragraph: 

In a practical biology laboratory, it is not necessary to maintain a dust-free 
environment. Hence, fine microscopic dust particles are nearly impossible to avoid. On 
laser excitation, these particles fluoresce to give high intensity spikes. Moreover, in some 
cases, bad mixtures of cDNA solutions result in precipitation, and these particles 
fluoresce with a very high intensity. These effects are simulated by adding spike noise at 
a preset rate. Such intensity spikes are added randomly across the entire slide area, the 
number of such noise pixels being preset in terms of the total number of pixels in the 
array. The amount of spike noise in an array is set with reference to the percentage, L^pi, 
of the total number of pixels in the array. Typical low to high noise levels are be set by 
selecting 0.1% to 10%. Once a pixel is selected for spike noise, the adjacent pixels have a 
higher probability of being affected. Thus, a random number, Wspi, of pixels are chosen in 
an arbitrary direction to be influenced by this noise. The intensity, Ns, of the spike noise 
is governed by an exponential distribution with mean ^^p,. In FIGURE 12, the 
exponential mean is fixed but the spike rate is increased through the parts of the figure. 
FIGURE 12 shows increased spike noise levels L^. FIGURE 12(a) shows a level of 
0.1%. FIGURE 12(b) shows a level of 5%. FIGURE 12(c) shows a level of 10%. An 
exponential rate range is maintained for the above. 

Please replace the paragraph beginning at page 17, line 24, with the following 
rewritten paragraph: 

Physical handling of the array slides can result in surface scratches. These 
typically result in low intensity levels. Scratch-noise intensity is parameterized as a ratio, 
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Ksc, giving the background-to-scratch-noise intensity level. Other parameters are the 
number of strips, strip thickness W^c, and a random strip length, Lsc, given as a multiple of 
the spot size. The latter is modeled as a uniform distribution: Lsc ~ U[Lsch Lsa]- Strips are 
placed at random positions on the array, and are inclined according to a (discrete) 
uniformly random angle, 9,^ g {0°, 45°, 90°, 135°, 180°}. FIGURE 13 shoves scratch 
noise with its parameter settings; the number of strips is fixed at 7. FIGURE 13 shows 
the noise for incremental parameter settings: (a)^^ ^ • ' U[2,J Jir ^ ~ 2.0, W^ ^ - 4 pixels; 
(b)-fc ^ ^5„10], ic^ -^=^^tOH^ ^ - 7 pixels; (c) ■ ■ U[1J51 ic^ -^4vOr^ ^ - 10 pixels. 
The number of strips is fixed at 7. FIGURE 13 (a) shows ~ L/F2 7], k^- 1.5, = 3 
pixels: FIGURE 13(b) shows L,. ~ U\5 151, - 2.5, - 7 pixels; FIGURE 13(c) 
shows Lv. ~ LT8 451, = 4.0, = 15 pixels. The number of scratches is maintained to 
7. The noise factor k,.. = 0.1. 

Please replace the paragraph beginning at page 18, line 5, with the following 
rewritten paragraph: 

Fine fabric dust particles on the slides can create snake-tailed strips on laser 

excitation. These strips are normally higher intensity than the signal level. To simulate 
this noise, an equiprobable multi-directional snake noise has been generated consisting of 
some number, N^eg, of segments. Analogously to scratch noise, the intensity 
parameterized as a ratio, k^„, giving the average-signal-to-snake-noise intensity level, the 
number of snakes, snake thickness Ws„, and a random length, Ls„, given as a multiple of 
the spot size. The latter is modeled as a uniform distribution: Ls„ ~ U[Ls„i, Ls„2]. FIGURE 
14 shows different parameter settings for snake noise. F IGURE 14 shows the noise for 
incremental parameter settings: FIGURE 15 (a) shows N ,.^ 5, Ls„ ~ U[5, 10], k,„ = 0.50, 
Ws„ = 2 pixels; FIGURE 15 (b) shows N^eg = 10, L,„ ~ U[5, 30], k,„ - 0.33, Pf,„ = 3 pixels; 
FIGURE 15 (c) shows N,eg = 1 5, L,„ ~ U[iS 5, 80], k,„ = 0.25, = 5 pixels. The 
direction of the tail was randomly chosen with equal probability for each. 
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Please replace the paragraph beginning at page 24, line 29, with the following 
rewritten paragraph: 

FIGURE 17 shows two microarrays generated with NS^ = 35 rows and NS^ = 25 
columns, at 5/, = 5^, 900 pixels per block. The example shows full size arrays 
simulation with different parameter settings. Depending on the parameters the arrays are 
called as "average" and "noisy" in quality. FIGURE 17(a) shows good quality, has SNR 
of 2.0, with normal background, and spike noise = 0.3%. Array boundaries are set at 
(Mt, Ml, Mr, Mt)= (100, 100, 100, 100). By choosing parameters, two different array 
qualities have been generated. Part a FIGURE 17(a) illustrates an ideal microarray image 
with normal background and parameters P = 3000, SNR = 2.0, a = 0.05, Gsp= 6, Pd= 
0.05, (da, db) = (2, 15), {ki^,h2) = (10, 10), Pou,iier= 0.05, L,^, = 0.3%, 5,^ =0.3, 

Please replace the paragraph beginning at page 25, line 13, with the following 
rewritten paragraph: 

FIGURE 17(b) shows a noisy array with SNR of 1.1 with parabolic background noise, 
and spike noise Lc ^ ; = 15%. FIGURE 17 (b) illustrates a noisy microarray image with 

parabolic background and parameters: (3 = 3000, SNR = 1 . 1 , a = 0.25, G,p= 4, Pq = OA, 
(da, db) = (15, 100), (ki^,ki,2) = (25,2 5), PoutUer^ 0.7, Lspi = 15%, 6ed = 0.03, 

Please replace the paragraph beginning at page 25, line 26, with the following 
rewritten paragraph: 

To illustrate how the simulation can be used to analyze microarray image 
software, the Array Suite [11] software is applied to extract the image intensities and 
ratios from the image and then these are compared to the corresponding intensities and 
ratios used for simulation. The ideal case is used to illustrate the utility of the simulation. 
FIGURE 18 shows a comparison between simulated signal (ideal setting) vs. extracted 
signal from a microarray image analysis program. In FIGURE i%& 18(a) , intensities from 
one fluorescent channel have been exfracted (y-axis) and plotted against the simulation 
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signal intensities. FIGURE 18(a), shows signal extracted from one fluorescent channel 
(y-axis) comparing to the signal used for simulation in the same channel (x-axis). T he 
extracted signal generally corresponds well to the simulated signal, with some variation. 
After excluding intensities less than 300, the mean and standard deviation of the 
difference between the two logio-transformed intensities are 0.016 (or 10° °'^ = 1.038) and 
0.038 (or 10° °^^ = 1.09), respectively. The ratio comparison is given in FIGURE 
■j-8b l8(b) . FIGURE 18(b) shows ratio from microarray image analysis program (y-axis) 
comparing to the ratios generated by the simulation (x-axis). When signal intensity is 
weak (less than 300), various noise components in the simulation process affect the 
accuracy of the signal extraction program. Since the problem may be unavoidable, a 
measurement quality metric may be used to provide confidence in downstream data 
analysis. In this case, if the signal intensity is less than 300, then the noise interaction is 
significant. 

Please replace the paragraph beginning at page 26, line 19, with the following 
rewritten paragraph: 

The simulation program has been used extensively in the design of the microarray 
image-analysis program used at the National Human Genome Research Institute. This has 
been done by testing the accuracy of the analysis program on simulated images exhibiting 
troublesome noise conditions and then tuning the program to achieve better results. One 
such application concerns large and overlapping spots, as illustrated in FIGURE 19(a), 
which shows part of an actual hybridized image in which some spots are substantially 
larger than intended owing to randomness in the cDNA deposition procedure. FIGURE 
19(a) shows part of an actual hybridized image with spots larger than average. This 
defect causes various problems, one being poor background estimation. This problem is 
illustrated by simulating an image with large spot size variation and drifting conditions 
[FIGURE 19(b)]. FIGURE 19(b) shows a simulated microarray with larger spots and 
spots overlapping with their neighbors. If the image analysis program extracts the local 
background by averaging the region around the bounding box (which was used as a 
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Starting condition in an earlier version of the NHGRI program), an elevated background 
average may be obtained since the bounding box may overlap neighboring targets that are 
large in size and strong in expression level. An additional problem is that some weak 
targets may not be detected [FIGURE 19(c)]. FIGURE 19(c) shows that an original 
background intensity extraction program produces an undetected spot (target in the 
middle without any outer boundary). Based on these considerations, the program has 
been modified to calculate the four average intensities from the four comers and the four 
average intensities from the four sides of the bounding box, and then take the minimum 
among these as the initial estimation of the local background. A histogram-based method 
is then invoked around the initial estimated background to fiirther improve the estimation 
[11]. The output from FIGURE 19(b) according to the modified program is shown in 
FIGURE 19(d): the weak target is detected and there is improved local background 
estimation for the spots. FIGURE 19(d^ shows that the improved background extraction 
program more accurately measures the local background intensity and effectively allows 
detection of weak targets. 

Please replace the paragraph beginning at page 28, line 29, with the following 
rewritten paragraph: 

13. See wmv.imgresearch.com . genome www. stanford.edu/microarray, 



J.G. Morton, et al., Three abundance classes in Hela cell messenger UNA. 1974, Nature, 
250(463): 199-240 

Please delete the paragraph beginning at page 29, line 3 1 to remove the reference 
to the Bishop reference. 



search. 



■nutecscience^ 



Bishop, J. O., 
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